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Section  1.  Introduction 

The  dynamics  of  a  large  class  of  engineering  systems  are  described  only 
implicitly,  for  instance,  those  of  non-linear  circuits, swing  dynamics  of 
an  interconnected  power  system,  as  also  thermodynamic  systems  far  from 
equilibrium.  The  implicit  definition  of  their  dynamics  is  as  follows: 
the  state  variables  are  constrained  to  satisfy  some  algebraic  equations, 
i.e.  they  are  constrained  to  lie  on  a  manifold  M  in  the  state  space. 

The  dynamics  on  this  manifold  M  are  then  specified  implicitly  by  specify¬ 
ing  only  the  projection  of  the  vector  field  on  M  onto  a  certain  base  space 
above  which  M  lies.  (i.e.  a  subspace  of  the  original  state  space  of  the 
same  dimension  as  M) .  The  process  of  obtaining  the  system  dynamics  explicitly 
consists  of  'lifting'  the  specified  velocities  onto  a  vector  field  on  M 
^lifting  is  the  inverse  of  projecting).  Lifting  nay  not,  however,  be 
possible  at  points  where  the  projection  map  (restricted  to  the  rancent 
space  of  the  constraint  manifold)  has  singularities.  This  singularity 
is  typically  resolved  by  regularization ,  i.e.  by  interpreting  the  algebraic 
constraint  equations  as  the  singularly  perturbed  limit  of  'parasitic'  or 
fast  dynamics.  The  dynamics  of  the  original  system  are  obtained  as  the 
degenerate  limit  of  the  dynamics  of  the  regularized  system  -  the  resulting 
trajectories  may  be  discontinuous  and  this  is  referred  to  as  jump  behavior. 

The  foregoing  deterministic  theory  needs  to  be  extended  to  a  stochastic 
theory  for  three  reasons: 

a)  The  conditions  under  which  the  limit  trajectories  to  the  regularizations 
exist  are  extremely  restrictive  so  as  to  exclude  several  systems  of 


interest . 
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b)  Frequently,  the  algebraic  constraint  equations  arise  from  the  macro¬ 
scopic  aggregation  of  microscopically  fluctuating  dynamics,  e.g.  the 
flow  of  current  in  a  resistor,  the  demand  for  electrical  power  at  a 
distribution  point  in  an  electrical  power  network.  More  generally, 
deterministic  equations  describing  thermodynamic  systems  are  of  this  kind. 
Thus,  the  algebraic  constraint  equations  contain  in  addition  a  rapidly 
fluctuating  (or  white  noise)  component. 

c)  The  methods  of  analysis  for  deterministic  systems  of  the  implicitly 
defined  kind  involve  techniques  of  bifurcation  theory  -  their  conclusions 
are  extremely  sensitive  to  imperfections  and  the  addition  of  white  noise. 

Since  in  all  the  situations  of  interest  to  us,  the  intensity  of  the 
additive  noise  is  small,  we  study  in  this  paper  the  dynamics  of  implicitly 
defined  dynamics  in  the  presence  of  small  additive  noise.  In  fact,  we 
compare  the  conclusions  of  the  stochastic  theory  with  those  of  the 
deterministic  theory  in  the  limit  that  the  noise  intensity  tends  to  zero. 

The  foregoing  process  requires  the  computation  of  two  sets  cf  limits: 
the  limit  that  the  regularization  tends  to  zero  and  the  limit  that  the 
intensity  of  the  additive  white  noise  tends  to  zero.  In  general,  these 
limits  do  not  commute.  We  explore  in  this  paper  the  modelling  issue  of 
which  sequence  of  limits  is  appropriate  in  the  context  of  a  specific  system. 
The  layout  of  the  paper  where  we  carry  out  this  program  is  as  follows: 

In  Section  2,  we  review  briefly  the  dynamics  of  deterministic 
constrained  systems  and  their  jump  behavior.  With  some  minor  modifications 
we  follow  here  our  earlier  work  [11]  and  the  references  contained  therein. 

In  Section  3,  we  begin  the  study  of  noisy  constrained  dynamical 
systems.  For  the  initial  study  we  use  as  tools  the  work  of  Papanicolaou, 
et  al  [1C]  on  martingale  approaches  to  limit  theorems.  To  study  the 
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dvnamics  of  noisy  constrained  systems  in  the  presence  of  sa*ll  noise, 
we  develop  and  use  in  our  context  Laplaces  method  of  steepest  descent. 

We  study  in  several  separate  cases,  the  comparison  between  tne  detersninist ; 
and  small  noise  theory,  describing:  (i)  now  the  stochastic  theory 

yields  conclusions  about  system  dynamics  when  tne  deterministic  theory 
fails  and  (ii)  how  the  pump  behavior  of  systems  is  modified  by  the 
presence  of  small  noise.  This  section  is  a  considerable  extension  of  out 
previous  work  in  the  context  of  phase  transitions  in  van  ici  VaaJs  gases 
[12],  Several  examples  are  presented  to  instantiate  our  results. 

In  Section  -l,  we  present  the  detailed  deterministic  analysis  of 
Section  2  applied  to  the  dynamics  of  an  emitter  coup  led  relaxat i or 
oscillator  circuit.  We  then  show  that  the  experimental  :w. -  lu*  ic.na  of 
Abidi  (11  on  the  dynamics  of  these  circuit!.  ;r.  tnc  presence  :t  -.mall  r.c  i :» 
see r  not  to  aoree  with  the  stochastic  theory  p  resented  ;r. 

In  Section  5,  we  discuss  the  sequences  of  limits  imtlied  by  t'.e 
development  of  Section  3  -  and  the  nature  cf  systems  f:r  which  tha-  Jcv<-.  >j 
ment  yields  the  correct  conclusions.  In  particular,  we  shew  that  the 
development  of  Section  3  is  relevant  to  systems  where  the  separation  ;r. 
time  scales  between  the  slow  and  fast  components  is  very  iarce  and  is 
more  important  than  the  small  intensity  of  the  white  noise  (characterized 
by  a  certain  sequence  of  limits  subsumed  by  Section  3)  -  for  instance  m 
phase  transitions,  reaction  rates  and  other  phenomena  of  oon-equi 1 itr lum 
thermodynamics.  For  non-linear  circuits,  however,  the  separation  m  time 
scales  is  less  marked,  so  that  we  present  here  the  relevant  analysis 
(sample  function  calculations)  for  these  systems  with  the  order  of  limits 
reversed  from  that  of  Section  3.  We  use  as  tools  the  foundational 
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Note  that  cot  all  * >0  and  x  S  H  n  the  critical  points  (with  respect  to  y 
or  p  (x,y)  ate  the  equilibrium  points  of  the  deterministic  system  12.4) 
with  x  frozen  given  m  this  instance  by 

•  } 

V  *  -  t  grad,  Slx.y)  .  (3.11) 

•  y 

Further,  if  for  some  x  ,  Six  ,y)  is  a  Horse  function  (of  y) ,  then  for  all 

o  o 

_! 

■  >  3  every  local  maximum  of  p  (x,y>  is  a  stable  equilibrium  of  (3.11). 

To  compare  tne  noisy  constrained  system  with  the  deterministic 
.'onstratned  system  in  the  limit  that  V.  0,  it  will  be  necessary  to 
evaluate  integrals  like  (3.9)  in  the  limit  that  '•  •  3.  This  is  done 
using  tr.c  following  version  of  Laplace's  method: 

Theorem  3..  i Lap l a  re ’ 3  Me  thod ' 

n  •  •  • 

Let  for  each  x«$»  ,  Six.y)  have  global  minima  at  y.  vX) ,  (x) , . . ,y  (x) , 

where  N  nay  ie:end  or.  x.  Let  them  all  be  nor.- degenerate .  Further, 
let  S  x . y '  have  at  least  siairatie  growth  .  ir.  y;  as  y  *  ’.  Then,  in 
the  limit  that  •  ',  r  (x,y)  converges  to 


•i  .  /  '1 

a  ix!  5  (y  -•,■'/  a,  (x! 
i'-  1  ‘  /  i=l  1 


where  a . t x)  =  det(D*  3(x,y  |xl ) )  1 ' 

l  -  » 


(3.12) 


More  precisely,  if  t ( x , y )  is  a  smooth  function  having  polynomial  growth 
as  y  -  T,  then 

r 

)  p  (x,y ) dy 


litn  3.  x) 
’•  -  3  ’ 


=  lira  I  1  ( x , v 

no  J *m 
•f  •  / - 

=  a.  (x)  t  (x,y  (x)  )  /  a  (x) 

1=1  /  1=1 


(3.13) 


=  i  (X) 


the  Hessian  Dt  3(x,v)  at  v  =  v  (x)  is  nonsinqular 

2  •  ’l 


J 


-15- 


Proof:  Since  p  (x,y)  =  exp 


we  will  first  evaluate 


s (x,v)  .  , 

exp  -  — — * L- dy)  , 


L 


>>  (x,y)  exp  -  -  dy 


for  simplicity  first  assume  that  S(x,y)  has  a  single  global  minimum  at  y. 
We  will  then  show  that 


L 


d(x, y)  exp 


-  S^*-^dy  =  £{x,y  )  (27rAm>/2)  exp  -  S(x,y  )/A[l  +  °(1)J 


2  *1/2 
[det  S(x,y 


i.  14) 


First,  by  the  Morse  Lemmatsee  for  e.g.  Milner  (16])there  exists 

'  _  m  . 

a  neighborhood  CJ  of  y  and  a  change  of  coordinates  ®  t  given  oy 

—  * 

y  =  v (v)  such  that  y  =  v{0)  and 


3  ( x ,  y )  =  S  ( x ,  y  ) 


T  -  ('/.)' 

2  i=l  1 


i3.15} 


Further,  outside  the  neighborhood  U  of  y,  S(x,y)  >  S(x,y  )  +  £ 


for  some  :  >  D  so  that 


/0(x,y)  exp  -  ^ - dy 

m 


exp  [  - /■*•--]  o ( ,\  *" )  (3.16; 

A 


-  S  ( x ,  y  ) 

»rv/  ^  '-*7  -  l 

]R  "'/U 

for  all  £  >  0.  Clearly,  then  (3.1o)  does  not  contribute  to  the  leading 
term  of  (3.14)  Consider  now 


f  6  (x,y)  exp  -  -S  dy 

u  *  r 

;p  _  Slx^ j 


exp  (  -  7  y2)  *  (X,Y  (y)  )  ,det  D  Y(y)’|dv 

m  i=l  1 


(3.17) 


Now,  standard  manipulations  with  Gaussian  distributions  yield  that 


-1 
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2\ 


_  _ 

L  Vi  >  *(y)dy 

i=l 


=  (2^\)m/2  (i<0) 


+  o  (A)  ) 


(3.1 8) 


Thus,  to  evaluate  (3.17)  we  only  need  compute  |det  D  y(0)i  •  Differentiating 
(3.15)  twice  with  respect  to  y  yields 


D2  S(x,y)  =  ((^)-1)T  (^)_1 

2  dy  dy 


From  (3.19)  it  follows  that 


Idet  D  y (0) ! 


2  *  -1/2 
fdet  D2  S ( x , y  ) )  7 


(3.19) 


so  that  (3.14)  now  is  immediate  on  combining  (3.16),  3.17)  and  (3.16) 

In  the  instance  that  Svx,v)  has  several  global  minima 

*  *  * 

'.x!  ,  y,(x),  .  .  .,y  (x)  it  follows  from  an  easy  extension  of  the  fore¬ 
going  argument  that 


(2-0  "/2 


exo 


-S(x,v  ) 


N 

(  7  det 

i=l 


*  -1  2 

S  '.x,y  (x) )  ] 


:■  (x,y*  (x))+  ■'(!)] 


(3.20) 


Selling  7(x,y)  =  1  in 

/-S (x , y)  , 

exp  - -f-  dy. 

m 


(3.20:  yields  the  corresponding  expression  for 

Combining  this  with  (3.20)  we  have  equation  (3.13). 

O 


Remarks :  (1)  If  the  growth  conditions  on  S(x,y)  and  ^  < x , y )  are  uniform  in 

x  for  x,  <  R  it  can  be  shown  that  for  p  >  1 

'f.  (x)  -  T  (x)  ’P  dx  -  0  as  ♦  0  .3  21) 

e  ■  ' 
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Proof :  Is  presented  in  Sastry-Hijab  [121 . 

Remark :  The  order  of  the  limits  is  peculiar  in  Theorem  (3.3).  If  the 

order  is  interchanged  i.e.  \  4-  o  first  and  then  Z  *  0  it  is  clear  that 
one  recovers  in  the  limit  the  deterministic  development  of  Section  2 
(with  the  minor  modification  that  x  has  an  additive  white  noise  terms. 

The  jump-behavior  or  the  y-variable  is  as  explained  in  that  section.  If, 
however,  c  4  0  first  and  the  >.  t  0,  the  jump-behavior  of  the  y-vanable 
is  somewhat  different,  as  we  now  elaborate: 

The  behavior  of  the  conditional  density  of  y  given  x  as  *  »  0  is  as 

in  Theorem  3.2:  the  y  variable  is  at  one  of  the  global  minima  of  S(x,») 

2  * 

with  probability  proportional  to  the  curvature  of  S(x,  )  ( (Det  C,  Slx,v  )) 

at  that  minimum.  Consider  first  the  case  when  the  minimum  is  unique. 

There  is  then  a  jump  in  the  y-variable  if  there  is  a  change  in  the  global 
minimum  of  S(x,*)  as  x  is  varied.  Points  of  jump  then  will  be  points  of 
appearance  and  disappearance  of  global  minima  of  S ( x , * ) .  This  is  m 
contrast  to  tne  deterministic  picture  of  Section  2,  where,  for  the  instance 
that  g(X,y>  is  of  the  form  of  (3.11),  stable  equilibrium  of  the  sped-up 
system  3  are  local  minima  of  Six,*)  and  points  of  bifurcation  are  points 
appearance  and  disappearance  of  local  minima  of  S(x,»). 

We  illustrate  this  with  an  example  -  the  van  der  Pol  oscillator  of 
(2.5),  (2.6)  with  added  noise.  Consider 

x  =  y  +  >U  g(t) 

Zy  =  -x  -  y3  +  y  +  *TF  (t) 

4  2 

y  v 

Here  S(x,y)  =  -xy  -  so  that,  in  the  limit  that  c  *  0:  the 

x-process  converges  to  one  satisfying 
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S  A*  .2  * 


=  f  V' 


where 


L_ 


5  .  /  5 

fo  X)  =  -  a.  x  f  x , y  X'  /  a. 

i=i  /  s=: 
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plan*  t f i  figure  10.  Sc**e  or  rise  features  of  this  .«,r'e  are  noted  bel ow: 


|  ,i  for  --iy*  '  '•'  *  *‘3B  tft*  ®^u4tior.  tl.’j  h»  three 

solutions.  while  for  '■  -*q*  and  '  v  Uq*  the  equation 

haa  only  one  solution. 


'll) 

'HD 


As  '■■  *  *  .  t  -  and  *#  V  »  -*,  i  »  3  as>«pt©t leal ly . 


.  V- 

7 he  values  "  -  „'t  *.  i  »  and  V«  -.j  a.  i  •  2 1. 

0  .  P  0  0 

are  the  points  if  bifurcati  n  equation  (4.') 


treated  as  the  bifurcation  pre>«t«r,  i.e.  at  thete  point*  it 


;»  n't  possible  to  i-lv#  (4.h)  for  i  as  a  function  of  V  locally 
and  or.  i  I'ae .  y .  7hcs«  p  tr.ts  nay  be  shewn  te  &e  points  of 


1  \i  b i*  ur 


.at.  r. . 


Returning  new  ti  the  full  lystea  -  '4.t>)  and  4.')  we  see  that 

■ ;  ntmur  us  i  lut:->n,n  f  r  •>«  »y  stcr  ex;.*t  so  1  ■  no  as  i  car.  be  solved 
-  ntmu.  u..  ly  i  .  i  fur. -f  s  -n.  :  v  .n.  4.';  as  t  bt  air.: 


(4.0 


(4.6) 


When.  -:p  -  V„  •  ::n  -  ;>i  -  0 ,  i.e.  i  *  2R  or  i  *  21  ^  ~  vT.’4n  it 

appears  that  ^  as  infinite  sc  as  to  prevent  the  integration  cf  equations 
( 4 .  *i  1  ,  (4.4).  The  regularization  af  this  system  is  accomplished  by  taking 

int-  account  the  fact  that  parasitic  capacitances  present  in  the  transistors, 
is  well  as  the  finite  rate  cf  the  operational  amplifiers  will  prevent 

i  fr-m  varyir.ej  di  scent  mucus  ly  and  in  effect  change  the  description  of  the 
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of  the  circuit  dynamics  from  (4.6),  (4.7)  to 


dV  <VV> 


-  -  v  -  (2X0  -  21)  R  -  VTl:  (2I0-i)/l 


Equations  (4.6)  and  (4.9)  are  a  gross  simplification  of  all  the  actual 
parasitics  present  in  the  circuit.  A  more  detailed  and  exhaustive  description 
involving  all  the  parasitics  would  start  from  the  original  equations  (4.1)  - 
(4.5).  The  present  regularized  model  is,  however,  accurate  enough  for  our 
purposes.  The  phase  portrait  of  this  system  shown  in  Figure  1  includes 
a  single  unstable  equilibrium  point  (V=0,  i =1 ^ '  and  a  limit  cycle.  The 
limit  trajectories  of  (4.6),  (4.9)  as  :V0  exist  and  include  the  relaxation 
oscillation  shewn  in  Figure  12  -  l  limit  cycle  with  two  discontinuities  - 
at  the  points  where  the  tra-ectory  switches  from  the  ~1  on,  52  =  ff  'state* 
to  the  off,  Q2  or.  'state'  and  vice  versa.  Note  also  from  Figure  11 
chat  the  5*  on,  7-  cn  'state'  is  unstable  as  evidenced  by  the  trajectories 
of  (4.6),  (4.9)  ooir.tina  away  from  that  'state'.  The  current  waveform 
i(t)  is  as  shewn  ir.  Figure  13.  The  half  period  of  the  oscillation  T  may  be 
estimated  approximately  by  integrating  equation  4.8  with  the  approximation 


that  for  0£t_<T,i  «  1^,  so  that  we  have 


vt/:r 


(-2R  +  — )  di 

i 


T  =  ~~  [2R(-Vt/2R  +  Ij)  +  V?  .-MV^IjR)  j 


(4  JO  ) 
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From  equation  (4.10)  it  follows  that  the  frequency  of  oscillation  is 
(approximately)  lienarly  proportional  to  IQ,  which  enables  this  oscillation 
to  be  used  as  an  electronically  tunable  oscillator  (e.g.  in  a  phase  locked 
loop).  In  such  applications,  it  is  important  to  know  the  noise  characteristics 
of  the  oscillator  in  response  to  resistive  thermal  noise.  Experimental 
observations  of  Abidi  [1]  indicate  that  the  actual  (noisy)  cur  -ent  waveform 
is  as  shown  in  Figure  14.  Key  features  of  this  figure  are  as  follows: 

(a)  the  transitions  or  jumps  appear  to  be  noise  free 

(b)  the  noise  superimposed  on  the  deterministic  waveform  of 
Figure  13  appears  to  be  small  (low  intensity)  immediately 
following  a  jump  and  then  appear  to  build  in  intensity. 

We  assume  (see  •;.?.[14])  that  all  the  noise  sources  in  the  circuit 

can  be  lumped  into  a  single-r.oisv  current  source  i  (t)  shown  dotted  in 

n 

Figure  9:  i  (t)  is  assumed  to  be  white  with  intensity  X (with  X  small 
at  room  temperatures,  since  it  is  proportional  to  KT) .  It  is  easy  to  check 
that  the  equation  (4.6)  is  now  unchanged,  while  (4.7)  changes  to 

0  =  V  -  (21.  -  2i)R  -  V  Cr.(2i  -i)/i  +  2RvTT  (t)  (4.11) 

u  l  U  n 

V.'e  regularize  the  system  (4.6)  ,  (4.11)  as  before  to  obtain 

<  —  (Iq  -  i) /C  (4.6) 

£i  =  V  -  (21  -  2  i )  R  -  Vm  in  (21  -  i)  /i  +  2R  veA"  i  (t)  (4.12) 

Note  that  C  scales  the  intensity  of  the  white  noise  in  (4.12)  precisely 
for  the  same  reason  as  in  equation  (3.2)  of  Section  3.  The  techniques  of 
Section  3.2  may  now  be  used  to  obtain  that  as  c40,  the  V-process  converges 
weakly  on  C  (  f  0  ,T ] ;  IR)  to  one  satisfying  : 
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V  =  (I  -  i  A  (V)  )/C 

-  \ 

where  i  (v)  is  i  integrated  over  the  conditional  density  for  1  given 
V,in  the  limit  that  e+0,  p  \i  ,v).  As  in  the  example  of  Section  3.2,  we 
have  in  the  limit  that  XiO,  p\i,v)  converging  to  a  sequence  of  delta 
functions  jumping  from  one  leg  of  the  solution  curve  to  (4.7)  to  the  other 
atv=0.  Also,  choosing  the  interval  of  weak  convergence  to  be  large  it 
appears  that  the  relaxation  oscillation  is  borken  up. 

This  analysis  is  contrary  to  the  experimental  evidence  of  Abidi  [1] 
What  has  gone  wrong?  How  does  one  recover  the  experimental  results  of 
Abidi  [1]?  These  are  the  questions  that  we  taken  up  next. 
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Section  5.  Sample  Function  Calculations. 

The  mathematical  reason  for  the  anomaly  between  the  machinery  developed 
in  Section  3  and  the  experimental  conclusions  of  Section  4,  is  the  order  of 
limits  £40  followed  by  X40  in  Theorem  (3.3).  This  order  of  taking  limits  is 
suitable  for  explaining  phenomena  in  several  situations  in  non-equilibrium 
thermodynamics  (for  e.g.  phase  transitions  of  the  kind  discussed  in  Sastry- 
Hijab  [12],  Eyring  chemical  reaction  rates,  etc.  -  see  for  e.g.  Nicolis- 
Prigogine  [9],  Landauer  [6]).  In  fact,  it  has  been  noted  by  thermodynamicists 
of  the  Brussels  School  that "fluctuations  play  a  crucial  role  in  changing 
the  behavior  of  systems  near  bifurcation  fronts".  However,  this  order  of 
limits  is  not  fully  satisfactory  in  the  circuit  context.  The  reason  for 
this  lies  in  the  fact  that  the  order  of  limits  £-10  followed  by  X40  (Theorem 
3.3)  yields  the  correct  conclusions  only  when  the  dynamics  of  the  fast  (sped- 
up)  system  are  much  faster  than  those  of  the  slower  x  variable.  This  is 
so,  because,  as  we  state  in  Section  3.4,  Laplace's  method  of  steepest 
descent  picks  for  the  limit  values  of  F‘'"(x,y)  as  ’.*0  the  most  stable 
■-limit  sets  of  the  underlying  deterministic  systems.  This  in  turn  is 
consistent  with  the  intuition  that  in  the  presence  of  persistent  random 
perturbation  (wide-band  in  nature)  the  trajectories  of  a  system  will 
concentrate  after  sufficiently  long  periods  of  time  in  the  vicinity  of 
the  most-stable  sets.  However,  the  sufficiently  long  periods  of  time 
may  be  very  large  indeed.  It  is  possible  to  show,  for  example  in  the 

gradient  case  of  Section  3.2  that  the  average  time  required  to  excape  from 

k/\ 

a  stable  equilibrium  is  of  the  order  of  e  for  some  k>0  (see  for  eg. 


Schuss  [15],  Ventsel-Freidlin  (15]). 
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By  tikmq  limits  in  the  sequence  -;iO  follower  by  '■»0.  the  mp licet jcn 

-k/1  -k;  '• 

is  that  •  is  smaller  than  e  ,  i.e.  is  at  least  le  ),  sc  that  the 
fast  system  has  sufficiently  much  time  tc  concentrate  ;r.  the  vicinity  c# 
its  .-limit  sets.  This  is  frequently  the  situation  ;r.  r.cr.-equi  1 11  :  .  urn 
theromodynamics  where  the  slower  dynamics  are  frequently  assumed  tc  be 
'quasi-static'.  In  the  circuit  context,  however,  the  separation  of  time 
scales  between  the  slow  and  fast  variable  is  net  as  large  as  is  implied  by 
the  theorem. 

As  noted  in  the  remark  following  hhecrem  3.’;  if  the  order  of  limits 
is  interchanged  (i.e.  '40  and  then  SO),  or.e  recovers  the  deterministic 
development  of  Section  2.  Before,  we  further  elaborate  and  make  prcciso 
the  statements  of  the  previous  paragraph  we  indicate  hew  -r.e  analyses  sample 
functions  of  the  process  generated  by  '  ' ,  (2.2)  : r.  the  limit  that 

*10  followed  by  -SO.  The  ma-cr  tool  ••  -his  development  is  the  work  -f 
Ve.ntsel-rreidlin  [13]. 

We  consider  here  sample  functions  of  the  process  generated  by 


x  =  ?;:<,y)  +  *  *'  ",  x(0)  = 

'y  =  f(x,y)  +  -  -  v  (0)  =  y0  . 

with  precisely  the  same  assumptions  as  in  Section  3.  Let  .  *  i 

[0,T]  -*■  ipn  x  3Pm  be  a  map  from  the  interval  |0,T]  to  the  >:  ,  y 

with  ,■  (0)  =  x  ,  1  (0)  =  y  .  Define,  for  this  trajectory,  the 
x  u  y  0 

I  _  (  V)  by 


(5.1) 

(5.2) 


x  y 
space 

uncticnal 


I.  (.) 


(t) 


dt 


)  ) 


(5.  3) 
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Figure  5:  Visualisation  of  the  Trajectories  on  the 

Upper,  Center  and  Bottom  Sheets  of  the  Cusp 
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Figure  9:  Showing  the  Hopf  Befurcation  for  the  Example  3.6 
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2I0 


F i Turo  13; 


.-'bowing  th<?  Relaxation  Oscillation  in  the 
Emitter  Coupled  Oscillator 


